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ABSTRACT 

The first cosmological results from the ESSENCE supernova survey (|Wood-Vasev et al.ll200l are 
extended to a wider range of cosmological models including dynamical dark energy and non-standard 
cosmological models. We fold in a greater number o f external data sets such as the recent Higher-;? 
release of high-redshift supernovae (|Riess et al.ll2007l ) as well as several complementary cosmological 
probes. Model comparison statistics such as the Bayesian and Akaike information criteria are applied 
to gauge the worth of models. These statistics favor models that give a good fit with fewer parameters. 

Based on this analysis, the preferred cosmological model is the flat cosmological constant model, 
where the expansion history of the universe can be adequately described with only one free parameter 
describing the energy content of the universe. Among the more exotic models that provide good fits to 
the data, we note a preference for models whose best-fit parameters reduce them to the cosmological 
constant model. 

Subject headings: cosmology: observations — supernovae : general 



1. INTRODUCTION 

Ever since Type la supernova (SN la) measurements 
first indicated an accelerating expansion o f the universe 
(jRiess et all 119981: iPerlmutter et al.lll999f) . the focus of 
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cosmology has shifted dramatically. Over the last sev- 
eral years, the primary aim of many cosmological obser- 
vations has been to discover the reason for this acceler- 
ated expansion. The name we often give to the unknown 
cause of the acceleration is "dark energy." In this con- 
text, dark energy represents not only the possibility of a 
hitherto undiscovered component of the energy density 
of the universe but also the possibility that the standard 
models of gravity and/or particle physics require revi- 
sion. 

The ESSENCE ("Equation of State: SupErNovae 
trace Cosmic Expansion") supernova survey is an 
ongoing project that aims to measure the equation- 
of-stat e parameter of dark e nergy to better than 
10% dKrisciurias etalj 120051 : ISollerman etaLI 120061 : 
iMiknaitis et al.ll2007l ). In this paper, we use the first 
cosmological results paper from the ESSENCE super- 
nova survey (jWood-Vasev et aLll2007l hereafter WV07), 
where distances and redshifts for a large sample of newly 
discovered high-redshift SNe la are reported. Moreover, 
WV07 performed consistent light-curve fitting of not 
only the ESSEN CE sample of supernovae but also 
the local sampl e (jHamuv et alj 1996: Ri ess et all 119991 : 
Uha et all 120061 ) as well as the SN la d ata released 
by th e SuperNova Legacy Survey (SNLS; lAstier et al.l 
2006). Combining this with constraints from baryo n 
acoustic oscillations (BAOs; Eiscnstci n et al.1 [2005), 
WV07 placed constraints on the dark energy equation- 
of-state parameter of w — —1.07 ± 0.09 (statistical 
1ct)±0.13 (systematic) and on the matter density of 
Q m = 0. 267^0^18 (statistical la) for a flat universe. 

To extend the analysis presented by WV07 we have 

(1) added the 30 SNe la detected at 0.216 < z < 1.755 
by the Hu b ble Sp ace Telescope (HST) as reported by 
iRiess et al.l (|2007f l: 

(2) included further constraints from a wider range of 
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complementary observations; 

(3) allowed for more complex cosmological models by 
both relaxing the assumption of a flat universe and test- 
ing non-standard models inspired by new fundamental 
physics; and 

(4) applied model-comparison statistics to decide on 
the model best preferred by the current data. 

Today's SN la results are accom modated in what ha s 
become the concordance cosmology (|Spergel et alj [2003). 
flanked by constraints on the matter density, f2 m , from 
large-scale structure measurements and on the flatness of 
space from cosmic microwave b ackground (CMB) mea- 
surements (|Spergel et al.ll2007t ). The concordance cos- 
mology is dominated by dark energy, fl x w 0.7, with 
present evidence being consistent with Einstein's cos- 
mological constant, A (Einstein! Il917t IZel'dovicbl Il967t 
iPadmanabhanl [20031 ) . We refer to the model with a cos- 
mological constant as the standard model, or A model. 

However, uncertainty remains over whether the cos- 
mological constant is indeed the complete description 
of dark energy or whether the dark energy might have 
more complex behavior. This question is motivated by 
the enormous discrepancy between the theoretical pre- 
diction fc^Uiecosjnological constant and its measured 
value ( Weinberg] 1 19891 ) . That we apparently live in an 
era when J7a and f2 m are almost equal, known as the 
"coincidence problem," also suggests that we may have 
an incomplete cosmological model. Thus, a variety of 
suggestions for new physics have emerged. 

Some of these suggestions take the form of variations 
of the equations of general relativi ty (e.g., the Dval i- 
Gabadadze-Porrati (DGP) model; iDvali et alj l2000f ). 
while others invoke more compl ex and evolving form s of 
dark energy (e.g., quintessence: ICaldwell et alJll998ft . A 
basic way to explore more complex models is to parame- 
terize the dark energy by an equation-of-state parameter 
w, relating the dark-energy pressure, p, to its density, p, 
viap = wpc 2 . (Hereafter we set c = 1.) This parameter 
may be time variable and characterizes how the energy 
density evolves with the scale factor, a: p cx a~ 3 ( 1+w \ 

This paper is organized as follows. In Sect. [2] we discuss 
information criteria in the context of cosmological model 
selection. Section [3] details the data sets used, their indi- 
vidual systematics and assumptions, and our method of 
combining them. Section |4] describes each model in turn 
and assesses which is preferred by the data, the results 
of which are discussed in Sect. O 



2. MODEL SELECTION VS FITTING 
PARAMETERS 

Parameter fitting and goodness-of-fit (GoF) 20 tests 
alone are not effective ways to decide between possible 
models. These statistical measures are based on the as- 
sumption that the underlying model is the correct one. 
The x 2 statistic can test the validity of a particular 
model, but comparing relative likelihoods based on the 
X 2 values of different models does not properly account 
for the structural differences between them. 

In other words, x 2 statistics are good at finding the 

20 Defined as GoF = T(v/2, X 2 /2)/F(u/2), where I>/2, X 2 / 2 ) 
is the incomplete gamma function and v is the number of degrees 
of freedom. It gives the probability of obtaining data that are a 
worse fit to the model, assuming that the model is correct. 



best parameters in a model but are insufficient for decid- 
ing whether the model itself is the best one. One might 
be tempted to prefer the model that gives the best fit to 
the data, defined as the lowest x 2 ■ However, this does 
not account for the relative complexity of the models. 
To give a blatant example, a lOth-order polynomial will 
always give an equal or better fit than a straight line to 
any data set, but this does not mean that any of the ex- 
tra eight coefficients have any significance. It just means 
that a model with more parameters will generally give an 
improved fit (always, if the simpler model is a subclass 
of the more complex one). 

Moreover, even though many cosmological models can 
be expressed in terms of a "to" that describes the dy- 
namical behavior of dark energy, the different functional 
parameterizations used by different models mean that 
they are not referring to the same thing. Consequently, 
the value of w that the data pr efer is integrally related 
to the model used in a fit (e.g. lZhao et~alTl2007l ). This 
difficulty makes it unwise to compare different models 
by simply considering likelihood contours or best-fit pa- 
rameters. For example, the constant, w, that appears 
in the standard dark-energy model (Sect. I4.1.4|) is not 
the same parameter as the constant, wo, that appears 
in the w(a) parameterization of the variable dark-energy 
model (Sect. 14.2. ip . So if the best-fit value of wo drifts 
away from —1 it does not rule out w = — 1. 

Instead we turn to information criteria (IC) to 
assess the strength of models. These statistics favor 
models that give a good fit with fewer parame- 
ters. In this p aper we use t he Bayesian information 
criterion (BIC; Schwari 1978f) and th e Akaike infor- 
mation criterion (AIC; lAkaikd 119741 ) to select the 
best-fit models. ILiddld (|2004f ) examines the use of 
information criteria in the context of cosmological 
observations, and we follow his prescription here. 
Previous explorations into AIC and BIC in a cosmo - 
logical context include iGodlowski fc Szvdlowskil 12 005) , 



ISzvdlowski fc Godlowskil (120061)." ISzvdlowski et al 
(120061) . iMagueiio fc Sorkinl ([20071 ) . iMukheriee et al 
(120061 ) andlBiesiada (2001) 



The BIC (also kno wn as the Schwarz information cri- 
terion; lSchwardfl978h is given by 

BIC = -21n£ + fclniV, (1) 

where C is the maximum likelihood, fc is the number of 
parameters, and N is the number of data points used in 
the fit. Note that for Gaussian errors, x 2 — ~ 21n£, and 
the difference in BIC can be simplified to ABIC = Ax 2 + 
AfclniV. A difference in BIC (ABIC) of 2 is considered 
positive evidence against the model with the higher BIC, 
while a ABIC of 6 is considered strong evidence (jLiddld 

l200l . 

The AIC (IAkaikdll974l ). 



AIC = -21n£ + 2fc, 



(2) 



gives results similar to the BIC approach, although the 
AIC is more lenient on models with extra parameters for 
an y reaso n ably s ized data s et (In A" > 2). As mentioned 
in ILiddld (|2007l ). ISugiural (|1978l ) give a version of the 
AIC corrected for small sample sizes, AIC C = AIC + 
2fc( fc — 1)/(N — fc — 1), which is i mport ant when N/k < 
40 (|Burnham fc Andersonl 120021 12004 ). The correction 
is negligible in our case (~ 0.06). 
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Both tests can be applied to unrelated models (the 
simpler model need not be nested in the more complex 
model) . These criteria make an attempt to quantify Oc- 
cam's razor. When two models fit the data equally well, 
the simpler model (the one with fewer free parameters) 
is preferred. 

A poor information criterion result can arise in two 
ways: (1) when the model is a poor fit to the data, re- 
gardless of the number of free parameters; in this case 
a large reduced \ 2 value indicates that the model does 
not explain the data; and (2) when the data are too poor 
to constrain the extra parameters in the model; such a 
model would be disfavored by information criteria if a 
simpler model is available, but it may well be that with 
improved data the more complex model becomes pre- 
ferred. 21 Thus, information criteria alone can at most 
say that a more complex model is not necessary to ex- 
plain the data. In this paper, we find this to be the 
current situation for dynamical dark-energy models. 

The simple prescription for information criteria above 
is limited. A more in-depth analysis of the improve- 
ment gained by more complex models would not sim- 
ply count parameters, but would consider how much the 
allowed volume in data space increases due to the addi- 
tion of extra parameters (i.e. how much more flexible 
the model becomes) , as well as any correlations between 
the parameters. Bayesian model selection is a technique 
that takes this into account using Bayesian evidence, 
i.e., the average likelihood of a mo del over its prior pa- 
rameter ranges. ISaini et"afl (|2004l ) pioneered the use o f 
Bayesian evidence in cosmology and lLiddle et al.l (2006) 
analyze a variety of different parameterizations of dark 
energy, w, finding that the standard cosmological con- 
stant remains the favo r ed m odel. Other studies include 
lElgarpy & Multamakil (|2006l ). For this paper, we have 
used the first approximation provided by the IC without 
calculating the full Bayesian evidence. This simpler ver- 
sion is sufficient for our purposes, and when systematic 
errors dominate the uncertainty in the data (as has now 
become true for SN data sets), further statistical analy- 
sis becomes unwarranted. IC also require no assumptions 
for the prior or the metric on the space of model param- 
eters. 



~ 10%. With the recent release of the first 4 years of 
data and their cosmological implications (WV07), we are 
now in an excellent position to probe further cosmological 
models. 

Another large, high-quality supernova search, the 
SNLS, recently pu blished an extensiv e and homogeneous 
first-year data set (jAstier et alj [2006). In combining the 
two data sets we have leaned on the light-curve fitting 
performed by WV07, who fit all SNe from these diffcr- 
ent data sets wit h the same light-curve fitter, MCLS2k2 
(|Jha et alJ l2007t) . Following WV07, the uncertainty 
added due to the intrinsic diversity of SNe la is 0.10 mag 
(assuming a peculiar velocity uncertainty of 400 km s _1 ). 
We have used the data from Table 9 in WV07, including 
only those SNe for which the light-curve fits passed the 
quality cuts. That includes 60 ESSENCE supernovae, 57 
SNLS supernovae, and 45 nearby supernovae. 

We also incorporate the new data release of 30 SNe la, 
at 0.216 < z < 1.755 detected by HST a nd classified as 
"gold" supernovae by iRiess et ail ([2007) . Such high- z 
data are particularly useful for this analysis because the 
SN la constraints on the evolution of w are improved as 
the range of redshifts is extended. We adopted the local 
supernovae that these samples had in common in order 
to normalize the luminosity distances of the samples, and 
we included the uncertainty in this normalization in the 
distance errors for the HST SNe la. 

Ideally these two data sets should both be generated 
using the same light-curve fitter, in which case no nor- 
malization would be required. This is in progress, and in 
the interim we provide the combined data set as used in 
this paper. 22 

We note that any statistical analysis of the type we 
perform here may be limited by systematic errors in the 
data. Possible sources of systematic error in the super- 
nova data include local velocity structures (|Jha et al.l 



20071 IZehavi et al 



(jWood-Vasev et al 



1998) 
I2007P 



and the treatment of dust 



Astier et al.ll2006l ). WV0 7 con- 



centrate much of their discussion on the analysis of sys- 
tematic errors and how to minimize them. They calcu- 
late that the systematic error in supernova data gives a 
maximum uncertainty in w of 0.13 for the flat, constant- 
ly model when combined with BAO constraints. 



3. DATA 

In order to test the different models, we have used 
observational data from a variety of sources described 
below. 

3.1. Type la Supernovae (SNe la) 

As mentioned in the introduction, one of the prime 
new data sets in this work is the SN la data from the 
ESSENCE project. The ESSENCE project is a ground- 
based survey designed to detect and monitor about 200 
SNe la in the redshift range z = 0.2 — 0.8. 

The strate gy and implementatio n of this project are 
described bv lMiknaitis et all (|2007l ), and the goal of the 
completed survey is to constrain w with a precision of 

21 For example, general relativity would fail to rival Newtonian 
gravity if the only experiment available were dropping balls from 
the leaning tower of Pisa, but general relativity would become nec- 
essary when data on the precession of Mercury's orbit became avail- 
able. 



3.2. Cosmic Microwave Background (CMB) 

The characteristic angular scale of the first peak in the 
CMB anisotropy spectrum is given by, 



9a 



r s (zis) 
Da(zi s ) 



(3) 



where r s (zi s ) is the comoving size of the sound horizon 
at last scattering [roughly proportional to (Q m H^ ) -1 ' 2 ] 
and Da{zi s ) is the comoving angular distance to the last- 
scattering surface. 

Foll owing the prescr i ption given by iDoran fc Lillevl 
( 2002) and iPage et all (|2003l ). we have converted the 



The 



data 



available 



at 



http : //www. dark- cosmology . dk/archive/SN 

http : //braeburn.pha. jhu. edu/$\sim$ariess/R06 and 
http://www.ctio.noao.edu/essence When used, please cite 
WV07, iKiess et all [120071) and lAstier et all {2006) in addition 
to this paper. These Web sites will be updated when the 
self-consistent light-curve-fitting is complete. 
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WMAP three-year result (jSpergel et al.ll2007f) on the lo- 
cation of the first peak to a reduced distance to the last- 
scattering surface, 



R = 



in*: 



s k 



dz 
H(zj 



= 1.71 ±0.03, 
( 4 ) 

where Sk(x) = sina;,a;, and sinhx for k = +1,0, and — 1, 
respectively. 23 

The value of this parameter is somewhat model de- 
pendent. For example, it changes slightly when ma ssive 
neutrinos are included (|Elgar0v fe Mu ltamaki 2007]), and 
it is weakly dependent on the density of the dark energy 
at the last scattering surface (through the Ctf s term: see 
iDoran fc LilTevI l2002f ) . Using R = 1.71 therefore artifi- 
cially excludes some models with, for example, strongly 
varying dark energy. However, constraints from big bang 
nucleosynthesis (BBN) rule out models with vastly dif- 
ferent expansion dynamics from the standard model at 
the time of nuc l eosynthesis (|Carroll fc Kaplinghalj |2002t 
iSteigmanl [20061 : 1 Wright! [2001 7 so this is not an unrea- 
sonable exclusion. The BBN results may change if not 
only gravitational dynamics but also particle physics pro- 
cesses changed in any of the models. 

The robustness of t he sh ift parameter is tested in 
lElgarpy fc Multamaki (|2007| ) compared to fitting the full 
CMB power spectrum, and they confirm that it is an 
accurate measure for non-standard cosmologies such as 
those tested here. Degeneracies that arise from using R 
rather than fitting the full CMB power spectrum are well 
constrained by other data such as BAOs and SNe. 

3.3. Baryon Acoustic Oscillations (BAO) 

Similar to the use of the angular scale of the first peak 
in the CMB spectrum, we can use the measurement of 
the peak at ~ 100ft." 1 Mpc of the large-scale correlation 
function of luminous red galaxies in the Sloan Digital 
Sky Survey (jEisenstein et al.ll2005[ ) to further constrain 
cosmological parameters. 

The large-scale correlation function is a combination of 
the correlations measured in the radial (redshift space) 
and the transverse (angular space) direction, and thus 
the relevant distance measure is the so-called dilation 
scale, 



Dv(z) = [DAizfz/Hiz)} 1 ^ , 



(5) 



at the typical redshift of the galaxy sample, z = 0.35. 
The absolute scale of the BAO is given by the sound 
horizon at last scattering, and the dimensionless combi- 
nation A(z) = Dy {z)\J Cl m HQ I z is well constrained by 
the BAO data to be 



A(0.35) = Dy(0.35) 



0.35 



= 0.469 ±0.017. 



(6) 



23 Other 
IWang 



calculations of the 



shift parameter include 
Mukhcricc (200 3), who find R = 1.70 ±0.03, and this is 



used by, for exam ple^ Alam ct al. (20Q3) and ILiddle et al.l (2006) 
IElgar0v fc Multamakil l|2007l) find R = 1.71 ± 0.03. We calculated 
the v alue independently. In a new paper IWang & M ukherieel 
(2007) find R = 1.71 ± 0.03 when allowing nonzero cosmic 
curvature. The small difference has a negligible effect on the 
results. 



TABLE 1 

Summary of models 



Model 



Abbrev. a Parameters' 5 Section 



Flat cosmo. const. 


FA 




4.1.1 


Cosmological const. 


A 




4.1.2 


Flat constant w 


Fw 


Q m , W 


4.1.3 


Constant w 


w 




4.1.4 


Flat w(a) 


Fwa 


U m , WO, W a 


4.2.1 


DGP 


DGP 




4.3.1 


Flat DGP 


FDGP 


n ro 


4.3.2 


Cardassian 


Ca 


fl m , q, n 


4~4l 


Flat Gen. Chaplygin 


FGCh 


A, a 


4.5.1 


Gen. Chaplygin 


GCh 


Ofc, A, a 


4.5.1 


Flat Chaplygin 


FCh 


A 


4.5.2 


Chaplygin 


Ch 


n k , a 


4.5.2 



a The abbreviations used in Fig.[9] b The free parameters in 
each model. Note that when fitting the SN la data we also 
fit an additional parameter, M, for the normalization of SN 
magnitudes. We include this in the number of degrees of 
freedom and in k, but have not listed it here as a parameter 
in each model. 

3.4. Not Used 

There are additional sources of observations, which are 
potentially important for constraining cosmological pa- 
rameters, that we have not included in this analysis. We 
have, for example, chosen t o omit any information from 
distant gamma-ray bursts fGhirlanda et al. 2004) since 
such data are rather controvers ial ((Mo rtsell fc Sollerman! 
I200a iFriedman fc Bloomll2005|) . We have also not used 
the c onstraints from X -ray data on relaxed galaxy clus- 
ters ((Allen et al.ll2004T ). although this method may well 
become important in the future. The weak lensing 
data from the Canad a-France-Hawaii Telescope (CFHT; 
iHoekstra et al.l I2006T) are an additional source of data 
that we have omitted here. At least some of these ad- 
ditional sources of information will become increasingly 
important as the data improve. 

Apart from geometrical probes, as a consistency check 
we have used large-scale structure constraints on the 
growth factor as measured by the Two Degree Field 
Galaxy Redshift Survey. The constraints on the models 
from the growth factor are similar, but weak, compared 
to the BAO constraints. Since the computation of the 
growth factor in modified gravity models is potentially 
very complicated, we have chosen not to include those 
constraints in the final analysis. 

3.5. Combining the Constraints 

Both the distance to the last-scattering surface de- 
rived from the CMB and the BAO constraints depend 
on the baryon density (flbh 2 ) and its uncertainty, and 
for this we use the value obtained from WMAP. How- 
ever, the depen dence of BAO on bary on density is weak 
[oc (fl b h 2 )-°- Q8 (jEisenstein et al.ll2005l )]. and for distance- 
related parameters, BAOs and CMB are independent. 
Since the SNe la, CMB, and BAOs are effectively inde- 
pendent measurements, we can combine our results by 
simply multiplying the likelihood functions. 

4. MODELS 

Hitherto, all measurements of w have been consis- 
tent with a cos molog ical constant, w = —1 (e.g., 
Garnavich et all 119981: iTonrv et all l2003t iKnop et all 
2003tlHannestad fc Mor tsell 2004; Ei senstein et al.ll2005l : 
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Fig. 1. — Flat dark-energy model: a flat universe with constant w 
(Sect. [47TT3jl . The constraint from each of the observational probes 
is shown by shaded contours (according to the legend). These are 
all 95% confidence intervals for two parameters. The combined 
contours (95% and 99.9% confidence intervals) are overlayed in 
black. The complementarity of the different observational probes 
is clearly demonstrated in the differing angles of the overlapping 
contours. The combined data form a clear preference around the 
cosmological constant model (w = —1). Despite the extra freedom 
afforded by allowing the dark energy to have an equation-of-statc 
parameter that differs from —1, the data do not show any indication 
that this freedom is required. 



Astier et al.ll2006t ISpergel et all [20071 : iRiess et ail 120071; 
Wood-Vasev eTatl I2007D ." However, as discussed in the 
introduction, this is not an unproblematic conclusion. 
Given the theoretical difficulties in predicting a cos- 
mological constant with the right vacuum energy den- 
sity, a variety of suggestions for other new physics have 
emerged. Many models use evolving scalar fields: so- 
called quintessence models, which allow a time-varying 
equation of state to track the matter density. In such 
models, the time-averaged absolute value of w is likely 
to differ from unity. Many o f these have b e en tes ted by 
other authors; for example, IWilson et alj ((2006) show 
that one specific ve rsion of quintessence, the inverse 
power-law potential ([Peebles fc Ratralll988P ). has a best- 
fit solution consistent with the cosmological constant 
model. Many other models including all kinds of ex- 
otica have been proposed, such as fc-essence, domain 
walls, frustrated topologic al defects and extra dimen- 
sions. iPadmanabhanl ([20031 ) gives a review of dark energy 
and its alternatives. 

Models can be broadly classed into two groups: (1) 
those that invoke some form of extra component to the 
composition of the universe (such as dark energy or 
quintessence), or (2) those that invoke a variation in the 
equations governing gravity. In some cases the two are 
interchangeable descriptions of a single theory. 

Here we choose several of the most popular models 
discussed in the literature and examine whether they are 
consistent with the data currently available to us: 

1. Standard dark-energy models, including varying w. 

2. Dvali-Gabadadze-Porrati (DGP) brane world model. 

3. Cardassian expansion. 

4. Chaplygin gas. 

In the following sections we outline the basic equations 
governing the evolution of the expansion of the universe 
in each of the different models, calculate the best-fit val- 
ues of their parameters, and find their AAIC and ABIC 



values. The models used and the parameters that de- 
scribe each model are summarized in Table [T] For some 
relevant models we plot the likelihood contours of their 
parameters (see Figs.[TH6]). We show how the magnitude- 
rcdshift evolution of each model compares to the super- 
nova data in Figure [7] and how their best fits match the 
value of R measured from the CMB data and the value 
of A measured from the BAO data in Figure [8] The in- 
formation criteria results are summarized in Tableland 
shown in Figure [HI 

All the contours in Figures [T] - [6] represent two- 
parameter confidence intervals and the uncertainties 
quoted in the text are the 95% confidence level for one 
parameter. 

4.1. Dark-Energy Models with Constant Equation of 

State 

The dark-energy models with a constant equation-of- 
state parameter, w, are described by the following equa- 
tion relating Hubble's constant to the scale factor, a: 



H 2 



-.2 



n, 



(7) 



where VL m is the current value of the normalised matter 
density, the curvature of the universe is given by Q,k = 
1— fl x — firm an d &>x is the current value of the normalised 
dark-energy density. 

4.1.1. Flat, Cosmological Constant Model (Flat A) 

The standard cosmological model is the ACDM model, 
which invokes w = — 1 at all times, in the form of a 
cosmological constant, tt x — Q\. The simplest version 
of this model assumes a flat universe (CIa = 1 — TO ), so 

n a 

which only depends on one parameter. Our best-fit value 
is 

Q m = 0.27 ±0.04. 

This has the lowest AIC and BIC of all models tested, 
so AAIC and ABIC are measured with respect to this 
model (Tabled). 

4.1.2. The Cosmological Constant Model (A) 

Allowing for deviations from flatness allows one extra 
degree of freedom, 



Hi 



— — H + "A ■ 



a u a" 
(recalling Qf. = 1 — f2^ — ^m)- 



(9) 



4.1.3. The Flat Dark-Energy Model (Flat Constant w) 

Allowing instead for different values of w, but main 
taining the constraint of flatness, we have 



H 2 



a 3(l+™) 



(10) 



where ft x = 1 — f2 m . We plot the likelihood contours in 
Figure [TJ The best-fit parameters are, 



0.27 ±0.04, 



w 



-1.01 ±0.15. 
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Fig. 2.— Flat, variable w(a) model (Sect. HSU} . The contours 
are the same as in Fig. [TJ The parameters of this model are very 
poorly constrained by the current data. 

4.1.4. The Standard Dark-Energy Model (Constant w) 

Relaxing the constraint of flatness, we fit the most gen- 
eral form of constant-equation-of-state dark energy using 
all three independent parameters of Eq. [7J 

The AIC and BIC for these four models suggest that 
the simplest model adequately explains the data, and 
there is little evidence supporting the inclusion of extra 
parameters. 

4.2. Dark-Energy Models with Variable Equation of 

State 

Allowing the dark-energy equation-of-state parameter 
to vary as the universe evolves adds additional degrees 
of freedom to the model. The Friedmann equation for a 
varying w is given by Eq. [7] with the following replace- 
ment: 



o 3(i+«0 ^ exp 



1 l + w(a') 



da' 



(11) 



Below we discuss one possible parameterization of w, 
which is linear in scale factor. However, it is important 
to note that the time variation of w can be param- 
eteriz ed in many ways (e .g. . lHannestad fc Mortselll 
2001 iBarnes et al.1 pOOa iCalvo fc Marotol 120061: 



Huterer fc Peirid 120071 Riess et al.ll2007[ ). andchoosing 



a particular parameterization enforces a particular form 
of time evolution on the model, which may not be ap- 
propriate. Giving an analytic form to the time evolution 
of w can act like a strong prior and may give misleading 
results (for furth er discussion see iBassett et al.l 12004 
iRiess et al.l 120071 ). In what follows, we consider what 
has become the most-common parameterization of w 
but caution that a non-parametric approach could be 
preferred. Promising non-parametric techniques include 
a type of principal component analysis that allows the 
reconstruction of cosmological features such as p(z), 
w(z), and their derivatives, as a function of redshift 
("e.g.. iLinder fc Hutere^ 120051: IHuterer fe Cooravl 120051: 
iHuterer fc Peirisll2007t IRiess et al.ll2007[ ). 

4.2.1. Standard Parameterization (Flat w a ) 

Using the parameterization w(a) = w q + w a (l — a) 
(|Chevallier fc Polarskil 12001: iLinderl f2003h . Eq. 1111 sim- 



FlG. 3.— General DGP model (Sect. 14.3.10 . The dashed line 
shows the flat model. Here the contours from the different obser- 
vational constraints disagree and the model is thus strongly disfa- 
vored. 



plifies to 



3(l+tu ) 



a 3(l+w +w a ) e 3w a (l- 



(12) 

When we fit for this model we assume a flat universe, 
although loosening this constraint does not change the 
results considerably. The best-fit parameters are 



a m = 0.27 ±0.04, 



-o = -i.i±° ;t 



= 0.8 



+0.8 



The large flexibility of this model means that it is 
poorly constrained by current data (see Fig. [2]). 

This is also the model that is used in t he Dark Energy 
Task Force (DETF) figure of merit (FoM: lAlbrecht et all 
2006). The FoM is given by the inverse of the area of 
the 95% confidence interval in the wo-w a plane. The 
smaller the area (thus the larger the FoM) , the better the 
discriminating ability of the experiment. We calculate 
that the data sets used here have an FoM of ~ 1 . 

4.3. The DGP Models 



The DGP models (IDvali et ali kOOO) arise from a class 
of brane-related theories in which gravity leaks out into 
the bulk at large distances, resulting in the possibility of 
accelerated expansion. 

Notably, this theory provides for an accelerating uni- 
verse without adding any extra parameters, twro param- 
eters being sufficient to define the model. iLud (|2006f ) 
shows how the growth of large-scale structure proceeds 
in the DGP model, but the position of the BAO peak 
is not expected to be influenced by this modification 
(|Yamamoto et al.ll2006D . 

4.3.1. DGP Model 
The general DGP model is governed by the equation 




(13) 



where Q m = 1 — — 2a/TVJ-\/1 — flk- The parameter 
r c is the length scale beyond which gravity leaks out into 
the bulk, and fl r is related to this length scale by Sl r — 
l/(4r c 2 tf 2 ). 

As is evident from Figure O the overlapping region 
that is preferred by the supernova and BAO data seems 
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to be inconsistent with the CMB dat a. This concurs with 
the an alysis of the DGP model by iFairbairn fc Goobarl 
(2006). However, it should be noted that the model can- 
not be ruled out based on this observation alone since 
the GoF is ~ 0.2; that is, assuming that the underlying 
model is indeed DGP, the probability of finding an even 
worse fit is ~ 20 %. 



4.3.2. Flat DGP Model 

The flat DGP model can be considered a more con- 
strained version of the general DGP model. It has only 
one parameter to fit, fi ro , and serves as an illustrative 
example of the power of information criterion tests. The 
model gives a best-fit x 2 value of 210 for 192 degrees of 
freedom (dof). The GoF of this is 18%. The flat cos- 
mological constant model has a best-fit \ 2 value of 195 
for the same number of dof, giving a GoF of 44% (Ta- 
ble [2]). Thus, comparing the GoF for the models may 
not seem to warrant a strong preference for one model 
over the other. However, since the models have the same 
number of fitted parameters, the difference in the BIC is 
equal to the difference in the % 2 j ABIC » 14 (see Ta- 
ble [2] and Fig. [9]) , indicating a strong preference for the 
flat cosmological constant model. 

In order to assess the significance of the IC results, we 
have performed Monte Carlo tests with 1000 simulated 
data sets where the underlying cosmology is given by the 
flat DGP. We fitted flat DGP and flat A models to each 
simulated data set and then compared the BIC for the 
two fits. The ABIC obtained 24 is roughly Gaussian with 
ABIC = — 19±9. So if the underlying universe genuinely 
followed the DGP model, we would expect our measured 
ABIC to be negative. Instead we find a ABIC of +14. 
The highest value we obtain from our 1000 simulated 
DGP data sets is ABIC = 7, still far from the value of 
14 we obtain for the real data. When the underlying cos- 
mology is assumed to be flat A, our simulated data sets 
give ABIC = 18 ± 8, fully compatible with the measured 
value. From this we conclude that a large difference in 
the IC of different models indeed points to a very strong 
statistical preference for the model with the lower value 
of the IC. 

4.4. Cardassian Expansion 

Cardassian models (jFreese fe Lewis! I2002T ) involve a 
modification of the Friedmann equation that allows for 
acceleration in a flat, matter-dominated universe. The 
reason for the modification could be the self-interaction 
of dark matter, or the embedding of our observable 
three- dimensional brane in a higher-dimensional uni- 
verse. IWang et al.l (|2003f ) calculated that a Cardassian 
model could be discriminated from a generic quintessence 
model or a cosmological constant given expected future 
data se ts, such as the Sup erNova Acceleration Probe 
(SNAP; lAldering et aLT 2004) , one possible manifestation 
of the Joint Dark Energy Mission. We now employ the 
interim data sets that have become available to see what 
current data can determine. 

Note that we have assumed that any non-standard 
modifications to the location of the CMB and B AO peaks 
are negligible. 
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Fig. 4. — Cardassian expansion (Sect. l4.4j l. The dotted line shows 
the parameters that would agree with the flat, constant-?*) model 
(for a wide range of to-values). The cross marks the parameters 
that match the flat A model. This is a model with three free 
parameters (f2 m is not shown), and it is not very well constrained 
by the current data. 




0.80 



Fig. 5. — Flat generalized Chaplygin gas (Sect. I4.5.1D . The 
cross at a = 0, A = 0.72 marks the parameters that match the 
best-fit flat A model, while the dotted line shows the parameters 
that match the A model (with Q m = 1 — A). Again, despite the 
flexibility of this model, the best fit is achieved for parameters that 
are consistent with the flat A model. 



The original power-law Cardassian model has, 
H 2 Q m (1 — &m — &k) 



Hi 



~aJ 



-, 3vi 



(14) 



where n is a dimensionless parameter related to w. The 
original Cardassian model is equivalent to the standard 
dark-energy model (Sect. 14. 1 .4[) for w = n — 1, so there 
is no need to additionally fit that model. However, other 
incarnations of Cardassian expansion do not match any 
standard dark-energy model. One example is "modified 
polytropic Cardassian" expansion, which follows, 



H 2 



- - v ; 



1 



1 



a 3g(n-l) 



(15) 



ABIC = B1C H 



gJQflatDGP 



For q = 1, this collapses to the flat dark-energy model 
with w = ri — 1. 

Cardassian expansion fits the data well. This is due 
to its close phenomenological similarity with standard 
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0.40 



Fig. 6. — Standard Chaplygin gas (Sect 14.5. 2} . The dashed line 
shows the flat version of the model. Clearly this model is a very 
poor fit to the data. The subtleties of information criteria are not 
required to determine that this model is disfavored. 



dark-energy models. In particular, we note that the best- 
fit Cardassian expansion parameters are consistent with 
(within la of) those that make the Cardassian expansion 
collapse to one of the standard dark-energy models (see 
Fig.|4|). However, it suffers in AIC and BIC tests because 
of its larger number of parameters (three). 



4.5. Chaplygin gas 

Chaplygin gas models (Kamcnshc hik et al.l 1200 1[ ) in- 
voke a background fluid with p oc p~ a . They are moti- 
vated by brane-world scenarios (Bcnt o et al] |2002. and 
references therei n) and may be ab le to unify dark matter 
and dark energy (|Bilic et al.ll2002l ). We consider both the 
generalized and standard (a = I) Chaplygin gas models, 
with and without flatness constraints. 

4.5.1. Generalized Chaplygin Gas 

The generalized Chaplygin gas has an equation of state 
governed by p = —A/p a (with p > and A being a 
positive constant) and obeys 



H 



m a 2 1 



fi fc ) [A + 



a 3(l+a) 



(16) 



where the standard cosmological constant model is re- 
covered for a = and f2 m = (1 — Ofc)(l — A), The 
reduced distance to the last-scattering surface has been 
calculated as 



R 



' (i-n fc )(i- A) 
10*1 



Ho vl^j 



dz 



H(z) 



(17) 

No modifications for the location of the BAO peak have 
been made. The flat version requires Qk = 0. 

Out of all the non-standard cosmological models that 
we consider, Chaplygin gas models fare the best under 
the information criteria tests (see Table [2]), with the flat 
version slightly preferred (Fig. [5]). This is not unex- 
pected, as the best-fit pa rameters are again c lose to flat 
A (as also found in, e.g., iBean fc T5o re 2003). 




Fig. 7. — Best fit for each model, determined using supernova, 
BAO and CMB data, plotted against the supernova data. Distance 
modulus is shown with respect to the empty model. Here you can 
see that the models that provide good fits to the data (lines in 
the upper legend) are all essentially identical fits, while the models 
that could not mimic a cosmological constant fit poorly (lines in 
the lower legend). The grey points in the background are all the 
raw supernova data used in the fits, while for illustrative purposes 
only we show the binned data as l arge black fil l ed cir cles with two 
dimensional error bars. Following Riess et al. (2007) we use bins 
of nAz = 6, where n is the number of points in the bin and Az is 
the redshift range. Distance modulus error bars are the quadrature 
sum of the distance modulus uncertainties in the bin, while redshift 
error bars show the standard deviation of the redshifts in the bin. 
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Fig. 8. — Best-fit values of the CMB _R parameter and the BAO A 
parameter for each model. The solid lines and grey shaded regions 
show the measurement and lrr uncertainty for these parameters as 
measured using CMB or BAO. 

4.5.2. Standard Chaplygin Gas 

The standard Chaplygin gas (a = 1) has 



h 2 n k , 

w5 = -£ + fl 
H o a 



tit 



A 



(1-A) 



(18) 



Again we also test the flat version, which requires Qk = 0. 
These standard Chaplygin gas models may be the most 
basic models arising from d-brane theory, but they are 
not good fits to the data (Fig. see also IBean fc Dorel 
[2lM[Znul[200l . 

5. DISCUSSION AND SUMMARY 

We have tested a variety of non-standard cosmolog- 
ical models against the latest cosmological data. This 
includes new data of SNe la from the ESSENCE, SNLS, 
and Higher- z collaborations. We have also included the 
reduced distance to the last-scattering surface from the 
CMB and the constraints from BAO. Based on informa- 
tion criteria, both AIC and BIC, the simplest model of 
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TABLE 2 

Summary of the Information Criteria Results 





X I dof 


GoF (%) 


AAIC 


ABIC 


Flat cosmo. const. 


194.5 / 192 


43.7 








Flat Gen. Chaplygin 


193.9 / 191 


42.7 


1 


5 


Cosmological const. 


194.3 / 191 


42.0 


2 


5 


Flat constant w 


194.5 / 191 


41.7 


2 


5 


Flat w(a) 


193.8 / 190 


41.0 


3 


10 


Constant w 


193.9 / 190 


40.8 


3 


10 


Gen. Chaplygin 


193.9 / 190 


40.7 


3 


10 


Cardassian 


194.1 / 190 


40.4 


4 


10 


DGP 


207.4 / 191 


19.8 


15 


18 


Flat DGP 


210.1 / 192 


17.6 


16 


16 


Chaplygin 


220.4 / 191 


7.1 


28 


30 


Flat Chaplygin 


301.0 / 192 


0.0 


30 


30 



Note. — The flat cosmological constant (flat A) model is pre- 
ferred by both the AIC and the BIC. The AAIC and ABIC values 
for all other models in the table are then measured with respect 
to these lowest values. The goodness of fit (GoF) approximates 
the probability of finding a worse fit to the data. The models are 
given in order of increasing AAIC. 

a flat universe with a cosmological constant remains the 
best model to explain the current data. 

Information criteria provide a valuable way to get a 
relative ranking of the viability of scenarios, using a sta- 
tistical analysis that gives strong weight to the most sim- 
plistic model that fits the observations. This does not 
mean that the simplest model is always correct, rather 
that more complex and flexible models are not (yet) nec- 
essary. A poor information criteria result will arise when 
the data are not good enough to adequately constrain 
the model. In order to falsify a model one should look 
for contradictions in the data - such as when multiple 
data sets measure inconsistent values for the parameters 
of the model. This occurs, for example, in the standard 
Chaplygin gas model (Fig. [S]) . 

We provide a graphical representation of the IC results 
in Figured This shows not only the AAIC and ABIC, 
but also the number of parameters in each model (crosses 
and right-hand ordinate) . Given the current data, the flat 
cosmological constant model is clearly preferred by these 
tests. It almost achieves the best fit of all the models 
despite its economy of parameters. 

Following it are a series of models that give compara- 
bly good fits but have more free parameters. They are 
flat general Chaplygin gas, cosmological constant, and 
flat constant w, which all have two free parameters; and 
general Chaplygin gas, flat w(a), Cardassian expansion, 
and standard dark energy (constant w), which have three 
free parameters. We show how their magnitude-redshift 
evolution compares to the supernova data in Figure [7] 
(lines in the upper legend) and how well they fit the CMB 
and BAO data in Figure [8j Any of these models could 
eventually prove to be the best description of our uni- 
verse, but for the moment the data are not sharp enough 
to demonstrate the value of the extra complexity. Nev- 
ertheless, it seems suggestive that these models can all 
reduce to flat A and their best-fit parameters do so (to 
within la). The flat general Chaplygin model, for exam- 
ple, reduces to the flat cosmological constant model when 
a = and A = 1 — fi m . The actual values of the best fit 
are a = 0.03 ± 0.10 and A = 0.73 ± 0.04 (corresponding 
to n rn = 0.27 ±0.04). 
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Fig. 9. — Graphical representation of the results in Table [2] 
AAIC is represented by the light grey bars, ABIC by the dark grey. 
The order of models from left to right is the same as the order in 
Table [2] which is listed in order of increasing AAIC. The crosses 
mark the number of free parameters in each model (right-hand 
ordinate). The AIC and BIC qualitatively agree on the ranking 
of models, but the BIC is harsher than the AlC on models with 
more parameters. The primary increase in AIC for the first seven 
models (after flat A) is due to the larger number of parameters, 
whereas the last four models suffer because they are poor fits to 
the data. A ABIC of more than 2 (or 6) is cons idered positive (or 
strong) evidence against a model (Liddlc 2004). We express this 
as "unsupported" (or "strongly unsupported") by the data. 

Clearly new and better data are still needed to discrim- 
inate between these models. The advanced cosmological 
probes being pl anned for the next decade and beyond 
(jAlbrecht et afll200Gt IPeacock et al.ll2006D will improve 
considerably on current constraints and will be able to 
vigorously test the A model. If the results of these future 
experiments remain consistent with A, it would raise an 
interesting question. Is there a point at which we should 
accept A and abandon our scrutiny of more complex 
models? This is particularly problematic when alterna- 
tive models are able to match the A model to arbitrary 
precision. 

The answer is likely to be two-fold. Observers will 
continue to improve on current techniques and may dis- 
cover new techniques that could break the degeneracy for 
some of these models. Alternatively, theoretical consid- 
erations, such as the discovery of a quantum theory of 
gravity that makes accurate predictions in other realms 
of physics, may indicate a preference for a particular 
model. 

The last four models we tested, flat DGP, DGP, stan- 
dard Chaplygin, and flat standard Chaplygin, are clearly 
disfavored. They have fewer parameters than models like 
flat w(a), but they score poorly because they are unable 
to provide a good fit to the data. They do not reduce to 
flat A for any values of their parameters. 

In summary, given the current quality of the data, in- 
formation criteria indicate that there is no reason to pre- 
fer any more-complex model over the concordance cos- 
mology, the flat cosmological constant. It will be exciting 
to see whether future data sets change this conclusion. 
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